Evidence for an oscillatory signature in atmospheric neutrino oscillation 
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Muon neutrino disappearance probability as a function of neutrino flight length L over neutrino 
energy E was studied. A dip in the L/E distribution was observed in the data, as predicted from 



2 



the sinusoidal flavor transition probability of neutrino oscillation. The observed L/E distribution 
constrained i> M <-» v T neutrino oscillation parameters; 1.9 x 10 -3 < Am 2 < 3.0 x 10 _3 eV 2 and 
sin 2 29 > 0.90 at 90% confidence level. 

PACS numbers: 14.60.Pq, 96.40.TV 



Recent neutrino experiments using atmospheric [lj 



USUI, solar H0, 
accelerator neutrinos 



HElbj, reactor [R|, and 
have demonstrated that neu- 
trinos change flavor as they travel from the source to the 
detector, a phenomenon consistent with the hypothesis 
of neutrino oscillation. Neutrino oscillation is a natural 
consequence of neutrinos that have finite mass and fla- 
vor eigenstates that are superpositions of the mass eigen- 
states. The phenomenon is referred to as oscillation be- 
cause the survival probability of a given flavor, such as 



is given by: 



P{Vn -» z^) = 1 - sin 2 26 sin 2 



1.27Am 2 (eV 2 )L(km) 
E(GeV) 



(1) 

where E is the neutrino energy, L is the travel distance, 
Am 2 is the difference of the squared mass eigenvalues, 
and 9 is the mixing angle between flavor and mass states. 
This equation is true in vacuum for all cases, is true in 
matter for u T , but may be modified for oscillation 

involving v e which travel through matter. 

However, the sinusoidal L/E dependence of the sur- 
vival probability has not yet been observed. For solar 
neutrinos, the survival probability is non-sinusoidal as 
the two eigenstates in matter are no longer coherent after 
many oscillation cycles 15] . Reactor and accelerator neu- 
trino experiments have insufficient statistics at this time. 
The standard analysis Q, ^| of the large sample of at- 
mospheric neutrinos recorded by the Super-Kamiokande 
experiment has not been optimized to resolve the effect, 
although the zenith angle dependence strongly indicates 



maximal 



v T mixing with Am in the vicinity of 2 



to 2.5 x 1CP 3 eV 2 . The analysis described herein used 
a selected sample of these atmospheric neutrino events, 
those with good resolution in L/E. to search for the dip 
in oscillation probability expected when the argument of 
the second sine-squared term in Eq. \T\is ir/2. 

Super-Kamiokande (Super-K) is a cylindrical 50 kton 
water Cherenkov detector located at a depth of 2700 m 
water equivalent. The water tank is optically separated 
into two concentric cylindrical detector regions. The in- 
ner detector (ID) is instrumented with 11,146 inward fac- 
ing 20 inch diameter photomultiplier tubes (PMT). The 
outer detector (OD) is instrumented with 1,885 outward 
facing 8 inch PMTs. 

In the present analysis, 1489 live-day exposure of fully 
contained (FC) //-like and partially contained (PC) at- 
mospheric neutrino data were used. FC events deposit 
all of their Cherenkov light inside the ID, while PC events 
have an exiting particle that deposits visible energy in the 



OD. The direction and the momentum of charged parti- 
cles were reconstructed from the Cherenkov ring image. 
Each observed ring was identified as either e-like or /i-like 
based on the shape of the ring pattern. For FC multi-ring 
events, the particle type of the most energetic ring was 
used to identify /i-like events. Since more than 97% of 
PC events were estimated to be charged current (CC) 
interactions, all PC events were classified as /i-like. The 
atmospheric neutrino prediction in Super-K is modeled 
using a Monte Carlo (MC) simulation 

Event selection and classification in the present analy- 
sis are slightly different from those in the Super-K stan- 
dard oscillation analysis using zenith angle distributions. 
The fiducial volume for the FC events was expanded from 
22.5 kton to 26.4kton (event vertex should be more than 
1.5 m from the top and bottom walls of the ID and 1 m 
from the side wall) in order to increase the statistics of 
the data, especially of high energy muons. Estimated 
non-neutrino background in the expanded fiducial vol- 
ume was negligibly small, less than 0.1%. PC events 
were classified into two categories using the OD charge 
(photoelectron) information: "OD stopping events" and 
"OD through-going events" . Muons in the "OD stopping 
events" have stopped in the OD, while muons in the "OD 
through-going events" have passed through the OD. Fig- 
ure (1] shows the ratio of the observed charge in the OD 
to the expectation from the projected track length in the 
OD. Events with lower charge than the criterion were 
classified as "OD stopping events" . Since these two sam- 
ples have different resolution in L/E, different cuts were 
applied for each sample, improving the overall efficiency. 

The neutrino energy was estimated from the total en- 
ergy of charged particles observed in the ID. The en- 
ergy deposited in the OD was taken into account for PC 
events. The projected track length in the OD was used 
to estimate the energy deposited in the OD. The rela- 
tionship between the neutrino energy and the observed 
energy was determined based on the MC simulation. The 
flight length of neutrinos, which ranges from approxi- 
mately 15 km to 13,000 km depending on the zenith an- 
gle, was estimated from the reconstructed neutrino direc- 
tion. The neutrino direction was taken to be along the 
total momentum vector from all observed particles. The 
resolution of the reconstructed L/E was calculated at 
each point in the (cosO, E) plane, where O is the zenith 
angle. The energy resolution becomes poorer for higher 
energy PC events, due mainly to the saturation in the 
electronics that records the PMT charge. Therefore, ex- 
tremely high energy events (observed energy > 50 GcV) 
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FIG. 1: Observed charge (photoelectrons) in the OD divided 
by the expectation from the projected track length in the 
OD for the data (points), the OD through-going MC events 
(white region in histogram) and the OD stopping MC events 
(hatched region). The MC does not include oscillations and 
is normalized by the live-time. 
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Reconstructed zenith angle (cos©) 



Data MC + F M CC 



FC 


single-ring /i 


-like 


1619 


2105.8 


(98.3%) 




multi-ring /i- 


like 


502 


813.0 


(94.2%) 


PC 


OD stopping 




114 


137.0 


(95.4%) 




OD through- 


going 


491 


670.4 


(99.1%) 



TABLE I: Summary of atmospheric neutrino events used in 
the present analysis. Only /i-like events are used. Numbers 
of the MC events are normalized by the live-time. Neutrino 
oscillation is not included in the MC. Numbers in the paren- 
theses show the estimated fraction of i/ M + CC interactions 
in each sample. 



were excluded in this analysis. All the FC /z-like events 
have observed energy less than 25 GeV, so this cut is only 
relevant to PC events. Figure [5] shows 70% L/E reso- 
lution contours, as used for the selection criteria in this 
analysis. The reasons for the poor L/E resolutions are 
either large dL/dO for horizontal-going events or large 
scattering angles for low energy events. The bold solid 
central line in Fig. [3]} indicates the minimum survival 
probability of muon neutrinos predicted from neutrino 
oscillations with Am 2 = 2.4 x 1CT 3 eV 2 . It is clear that 
detecting high energy muon events is crucial to observe 
the hrst maximum oscillation in L/E. The resolution cut 
of A(L/E) < 70% was determined from the MC simula- 
tion to maximize the sensitivity to distinguish neutrino 
oscillation from other hypotheses. 

Table [I] summarizes the number of events used in this 
analysis after the L/E resolution cut. Figure shows the 
number of events as a function of L/E for the data and 
MC predictions. Two clusters of events are visible be- 
low and above 150 km/GeV. They mostly correspond to 



FIG. 2: Contour plots of 70% L/E resolution in the 
(cosO, E v ) plane for (a) FC single-ring, (b) FC multi-ring, 
(c) PC OD stopping and (d) PC OD through-going samples. 
Three additional lines in (a) show the survival probabilities 
of muon neutrinos predicted from neutrino oscillation with 
(sin 2 26>, Am 2 ) = (1.00, 2.4 x 10~ 3 eV 2 ). Full and half oscil- 
lation occur on the solid and dashed lines, respectively. 



downward-going and upward-going events, respectively. 
Below 150 km/GeV, the data and MC agree well. In 
Fig-Elthe data over non-oscillated MC ratio as a function 
of L/E is plotted together with the best-fit expectation 
for 2-flavor <-* v T oscillations with systematic errors. 
A dip, which should correspond to the first maximum 
oscillation, is observed around L/E = 500 km/GeV. We 
note that the position of the dip is about a factor of 3 
to 4 away from that of the predicted event number min- 
imum as seen in Fig. [31 Due to the L/E resolution of 
the detector, the second and higher maximum oscillation 
points should not be observable in this experiment. 

In order to confirm that the observed dip was not due 
to systematic effects, several tests were carried out. Sev- 
eral L/E distributions were made by changing the L/E 
resolution cut value. Plots based on the resolution cuts at 
60, 80 and 90% showed consistent dip structures as that 
based on the 70% cut. Also, L/E plots based on several 
other L/E bin sizes gave essentially the same results. In 
addition, the sign of the direction vector for each event 
was changed artificially. In this artificial data sample, 
the "upward-going" events should have little disappear- 
ance effect and therefore the L/E distribution should not 
show any dip structure around L/E = 500 km/GeV. The 
L/E distribution for this artificial data sample did not 
show any significant dip structure around 500 km/GeV. 
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FIG. 3: Number of events as a function of L/E for the data 
(points) and the atmospheric neutrino MC events without os- 
cillations (histogram). The MC is normalized by the detector 
live-time. 
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FIG. 4: Ratio of the data to the MC events without neutrino 
oscillation (points) as a function of the reconstructed L/E 
together with the best-fit expectation for 2-flavor «-» v T 
oscillations (solid line). The error bars are statistical only. 
Also shown are the best-fit expectation for neutrino decay 
(dashed line) and neutrino decoherence (dotted line). 



Finally, the L/E plot was made using FC single-ring e- 
like events. The e-like distribution was consistent with 
flat over the whole L/E range. Thus we are confident 
that the observed dip is not due to systematic effects in 
the event selection. 

The data/prediction at large L/E in Fig. 0] shows a 
slight rise from the expected flat distribution. We have 
studied possible causes of this deviation, and concluded 
that an energy-dependent systematic effects, such as the 
predicted neutrino interaction cross section, are the main 
sources of the non- flatness. The best-fit L/E distribu- 



tion for oscillations, allowing systematic terms to vary 
within the estimated uncertainty (as described below), 
also shows this rise with respect to no-oscillation predic- 
tion, as seen in the curves overlaid in Fig. ^ The rise at 
large L/E is consistent with the data. 

The observed L/E distribution was fit assuming «-> 
v r oscillations. The L/E distribution was divided into 
43 bins from log{L/E) = 0.0 to 4.3. The likelihood of 
the fit and the \ 2 were defined as: 

c{N ^ jN o bs) = g expHVf d )(ATf 



(2) 

(3) 
(4) 




7vr d =A?.PK^).(i+^. e 

3=1 

2 - _ 2ln f£( NPId > N ° hs )\ 



.))■ 



x 



\ C(N ohs , N ohs ) ) 



where N° hs is the number of the observed events in the 
i-th bin and iV l prd is the number of predicted events, in 
which neutrino oscillation and systematic uncertainties 
are considered. N® is the MC predicted number of events 
without oscillation for the i-th bin. Various systematic 
uncertainties are represented by 25 parameters tj , which 
include 7 uncertainty parameters from the flux calcula- 
tion (among these, absolute normalization is treated as 
a free parameter), 3 from the detector calibration and 
background, 2 from the data reduction, 5 from the event 
reconstruction, and 8 from the neutrino interaction sim- 
ulation. A more detailed description of the systematic 
error terms can be found in Ref. [lif. The second term 
in the likelihood definition represents the contributions 
from the systematic errors, where Uj is the estimated un- 
certainty in the parameter ej. The fractional effect of 
systematic error term Cj on the z-th bin is given by /j. 

A scan was carried out on a (sin 2 26, log Am 2 ) grid, 
minimizing y 2 by optimizing the systematic error param- 
eters at each point. The minimum % 2 was 37.9/40 DOF 
at (sin 2 26», Am 2 ) = (1.00, 2.4 x 10~ 3 eV 2 ). Including 
unphysical parameter region (sin 2 20 > 1), the best-fit 
was obtained at (sin 2 29, Am 2 ) = (1.02, 2.4 x 10~ 3 eV 2 ), 
in which the minimum \ 2 was 0.12 lower than that in 
the physical region. Figure El shows the contour plot of 
the allowed oscillation parameter regions. Three con- 
tours correspond to the 68%, 90% and 99% confidence 
level (C.L.) allowed regions, which are defined to be 
X 2 = Xmm+ 2 -48, 4.83, and 9.43, respectively, where 
Xmin i s the minimum \ 2 in the physical region. These in- 
tervals are derived based on a two dimensional extension 
of the method described in Ref. ,17]. The 90% C.L. al- 
lowed parameter region was obtained as 1.9 x 10 -3 eV 2 < 
Am 2 < 3.0 x 10~ 3 eV 2 and sin 2 29 > 0.90. The result 
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FIG. 5: 68, 90 and 99% C.L. allowed oscillation parameter re- 
gions for 2-flavor <-» v r oscillations obtained by the present 
analysis. 



is consistent with that of the oscillation analysis using 
zenith angle distributions [HH^. 

In order to test the significance of the dip in L/E, 
a null hypothesis that includes the basic shape of the 
L/E distribution is needed. The no-oscillation case 
was very strongly disfavored by the data at large L/E. 
We used two alternative hypotheses that basically re- 
produce the zenith angle dependent deficit, predict- 
ing that about half of muon neutrinos smoothly disap- 
pear at large L/E. These hypotheses are neutrino de- 
cay 0, E| and neutrino decoherence [2(3, EJ- The 
survival probability for neutrino decay is expressed 

as P(i/ p -» Vjj) = [sin 2 9 + cos 2 6»exp (-m/2r • L/E)] 2 
where r is the lifetime of a neutrino mass state. The neu- 
trino decoherence survival probability is Piy^ — > v^) — 
1— | sin 2 29 [1 — exp (—jqL/E)], where 70 is the decoher- 
ence parameter. Figure 0] includes the L/E distribution 
for the best-fit expectation for neutrino decay and neu- 
trino decoherence. The Xmin values were 49.1/40DOF at 
(cos 2 9, m/r) = (0.33, 1.26 x 10~ 2 GeV/km) for neutrino 
decay and 52.4/40DOF at (sin 2 29, 70) = 
1CP 21 GeV/km) for neutrino decoherence 
values were 11.3 (3.4 standard deviations) and 14.5 (3.8 
standard deviations) larger than Xmin f° r neu trino os- 
cillation. In order to check the statistical significance 
against the alternative models, 10,000 MC L/E distribu- 
tions were produced assuming neutrino decay with the 
best fit decay parameters. Each L/E distribution was 
fitted assuming neutrino decay and oscillation, and the 
X 2 difference for these two assumptions was calculated. 
For example, only 11 among 10,000 samples had x 2 for 
neutrino decay smaller by 11.3 or more than the same 



1.00, 1.23 x 
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rain 
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sample evaluated for neutrino oscillation. Therefore, the 
probability that neutrino decay could mimic neutrino os- 
cillations is approximately 0.1% as naively expected by 
3.4 standard deviations. The neutrino decoherence model 
is disfavored more strongly. 

In summary, we have studied the survival probabil- 
ity of muon neutrinos as a function of L/E using atmo- 
spheric neutrino events observed in Super-Kamiokande. 
A dip in the L/E distribution was observed around L/E 
= 500 km/GeV. This strongly constrains Am 2 . Alter- 
native models that could explain the zenith angle and 
energy dependent deficit of the atmospheric muon neu- 
trinos are disfavored, since they do not predict any dip 
in the L/E distribution. We conclude that the observed 
L/E distribution gives the first direct evidence that the 
neutrino survival probability obeys the sinusoidal func- 
tion as predicted by neutrino flavor oscillations. 
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